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This article describes the computation of cosmic microwave background anisotropies in a universe 

■ with multi-connected spatial sections and focuses on the implementation of the topology in standard 
' CMB computer codes. The key ingredient is the computation of the eigenmodes of the Laplacian 
\ with boundary conditions compatible with multi-connected space topology. The correlators of the 
. coefficients of the decomposition of the temperature fluctuation in spherical harmonics are computed 

^ ' and examples are given for spatially flat spaces and one family of spherical spaces, namely the 

fT^ ' lens spaces. Under the hypothesis of Gaussian initial conditions, these correlators encode all the 

fSJ , topological information of the CMB and suffice to simulate CMB maps. 

(N ■ 

CN ■ Preprint numbers: SPhT-T02/182, LPT-02/123, iastro-ph/0212223] 
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'. PACS numbers: 98.80.-q, 04.20.-q, 02.040.Pc 

o ■ 

O ■ I- INTRODUCTION 

^: 

^ ' Future cosmic microwave background (CMB) experiments such as the MAP (ij and later the Planck satellites 
• • , will provide full sky maps of CMB anisotropies (up to the galactic cut). These datasets offer the opportunity to 

■ probe the topological properties of our universe. A series of tests to detect the topology, including the use of the 
K> , angularpower spectrum jSj, |J, t5y6| , the distribution of matched patterns such as circles J} , correlation of antipodal 

j_] ■ points [3 and non Gaussianity |^ have been proposed (see Refs. 0, 0, 0| for reviews of CMB methods to search 
' for the topology) . The study of the detectability of the topology by any of these methods first requires simulating 
maps with the topological signature for a large set of topologies. These maps will allow one to test the detection 
methods, estimate their run time, and, once all sources of noise are added, determine to what extent a given method 
detects the top olog ical signal (in the same spirit as the investigation of the "crystallographic" methods based on 
galaxy catalogs |13|). 

In a simply connected^, spatially homogeneous and isotropic universe, the angular correlation function depends only 
on the angle between the two directions and the coefficients ag m of the decomposition of the temperature fluctuation 
in spherical harmonics, which are uncorrelated for different sets of £ and m. Multi-connectedness breaks the global 
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^ Geometers and cosmologists often refer to simple-connectedness as the "trivial topology". However, trivial topology has a different 
meaning in the context of point-set topology: in that formalism, the trivial topology is the smallest topology on a set X, namely the 
one in which the only open sets are the empty set and the entire set X. 
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isotropy and sometimes the global homogeneity of the universe, except in projective space (see, e.g., Ref. [T^l. 
Consequently, the CMB temperature angular correlation function will depend on the two directions of observation, 
not only on their relative angle, and possibly on the position of the observer. This induces correlations between the 
a^m of different t and m. Such correlations are hidden when one considers only the angular correlation function and 
its coefficients, the so-called Ci, in a Legendre polynomial decomposition, because they pick up only the isotropic part 
of the information and are therefore a poor indicator of the topology. This work aims to detail the whole computation 
of the correlation matrix 

^im =(a^ma^'m')i (1) 

which encodes all the topological properties of the CMB, and from which one can compute the usual Ci, simulate 
maps, and so on. 

The study of the detectability of a topological signal (if it exists) in forthcoming CMB datasets requires simulating 
high quality maps containing the topological signature for a wide class of topologies. Up to now, most CMB studies 
considered only compact Euclidean spaces [ESIsllEEiEl and some compact hyperbohc spaces [r^llMIT9ll20ll2lL[2^ . 
and focused mainly on the Ci. The approach developed in this article, and first introduced in Ref. jl4| . is well suited 
to simulate the required CMB maps in any topology once the eigenmodes of the Laplacian have been determined. It 
paves the way to the simulation of maps for a wide range of topologies, particularly spherical ones. 

Recent measurements of the density parameter $7 imply that the observable universe is "approximately flat"^, 
perhaps with a slight curvature. The exact constraint on the total density parameter obtained from CMB experiments 
depends on the priors used during the data analysis. For example with a prior on the nature of the initial conditions, 
the Hubble parameter and the age of the universe, recent analysis of the DASI, BOOMERanG, MAXIMA and 
DMR data [H HI IHl lead to n = 0.99 ± 0.12 at la level, and to f7 = 1.04 ± 0.05 at Icr level if one takes into 
account only the DASI, BOOMERanG and CBI data. Including stronger priors can indeed sharpen the bound. For 
instance, including information respectively on large scale structure and on supernovae data leads to = I.OIIqqq 
and n = 1.021°:°^ at 1(7 level while including both finally leads to = l.OOlgJg. This has been recently improved 
by the Archeops balloon experiments |2^ |2^ which get, with a prior on the Hubble constant, VL = 1.001qo2- 
conclusion, it is fair to assert that current cosmological observations set a reliable bound 0.9 < Vl < 1.1. These results 
are consistent with Friedmann-Lemaitre universe models with spherical, fiat or hyperbolic spatial sections. In the 
spherical and hyperbolic cases, ~ 1 implies that the curvature radius must be larger than the horizon radius. In 
all three cases — spherical, fiat, and hyperbolic — the universe may be simply connected or multiply connected. 

The possibility of detecting the topology of a nearly fiat universe was discussed in Ref. j^^. It was noted that 
the chances of detecting a multiply connected topology are worst in a large hyperbolic universe. The reason is 
that the typical translation distance between a cosmic source and its nearest topological image seems to be on the 
order of the curvature radius, and that when ~ 1 the distance to the last scattering surface is less than the half 
of that distance. See Refs. |28L l29l ISOl ISlL l33 | for some studies on detectability of nearly fiat hyperbolic universes. 
In a multiply connected flat universe the topology scale is completely independent of the horizon radius, because 
Euclidean geometry — unlike spherical and hyperbolic geometry — has no preferred scale and admits similarities. 
Note that in the Euclidean case, there are only ten compact topologies, which reduces and simplifies the analysis (in 
particular regarding the computation of the eigenmodes of the Laplacian) . In a spherical universe the topology scale 
depends on the curvature radius, but, in contrast to the hyperbolic case, as the topology of a spherical manifold 
gets more complicated, the typical distance between two images of a single cosmic source decreases. No matter how 
close 51 is to 1, only a finite number of spherical topologies are excluded from detection. The particular case of the 
detectability of lens spaces was studied in Ref. [s^l, which also considers the detectability of hyperbolic topologies. 

At present, the status of the constraints on the topology of the universe is still very preliminary. Regarding locally 
Euclidean spaces, it was shown on the basis of the COBE data that in the case of a vanishing cosmological constant the 
size of the fundamental domain of a 3-torus has to be larger than L > 4800 Mpc where the length L is 

related to the smallest wavenumber 2t:/L of the fundamental domain, which induces a suppression of fluctuations on 
scales beyond the size L of the fundamental domain. This constraint does not exclude a toroidal universe since there 
can be up to eight copies of the fundamental cell within our horizon. This constraint relies mainly on the fact that the 



^ The popular expression "flat universe" is misleading, because in General Relativity the "universe" is not three-, but four-dimensional 
and the Friedmann-Lemaitre solutions are dynamical, so that the universe is not flat — only its spatial section may be flat or nearly so. 
In what follows, we will implicitly assume we are talking about the three-dimensional spacelike sections of the universe when talking 
about flat, hyperbolic, or spherical spaces/universes. 
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smallest wavenumber is 27r/L, which induces a suppression of fluctuations on scales beyond the size of the fundamental 
domain. This result was generalized to all Euclidean manifolds in Ref. |l5j . A non- vanishing cosmological constant 
induces more power on large scales, via the integrated Sachs- Wolfe effect. For instance if JIa — 0.9 and fl^ = 0.1, 
the constraint is relaxed to allow for 49 copies of the fundamental cell within our horizon. The constraint is also 
milder in the case of compact hyperbolic manifolds and it was shown |l7l Hoj that the angular power spectrum was 
consistent with the COBE data on multipoles ranging from 2 to 20 for the Weeks and Thurston manifolds. Another 
approach, based on the method of images, was developed in pol[2ll l2 ^ . Only one spherical space using this method 
of images was considered in the literature, namely projective space jSjf. The tools developed in this article, as well 
as in our preceding works 0, IssL Is^ , will let us fill the gap regarding the simulation of CMB maps in spherical 
universes. 

Technically, in standard relativistic cosmology, the universe is described by a Friedmann-Lemaitre spacetime with 
locally isotropic and homogeneous spatial sections. These spatial sections can be defined as constant density or time 
hypersurfaces. With such a splitting, the equations of evolution of the cosmological perturbations, which give birth 
to the large scale structures of the universe, reduce to a set of coupled differential equations involving the Laplacian. 
This system is conveniently solved in Fourier space. In the case of a multiply connected universe, we visualize space as 
the quotient X/T of a simply connected space X (which is just a 3-sphere S^, a Euclidean space R'^, or a hyperbolic 
space H'^, depending on the curvature) by a group F of symmetries of X that is discrete and fixed point free. The 
group F is called the holonomy group. To solve the evolution equations we must first determine the eigenmodes TJ. 
and eigenvalues — of the Laplacian on X/T through the generalized Helmholtz equation 

ATifl=-4T[^l, (2) 

with 

4 = k^-K, (3) 

where k indexes the set of eigenmodes, the constant K is positive, zero or negative according to whether the space 
is spherical, flat or hyperbolic, respectively,^ and the boundary conditions are compatible with the given topology. 
The Laplacian in Eq. (0) is defined as A = D^Di, Di being the covariant derivative associated with the metric 
of the spatial sections {i, j = 1, 2, 3). The eigenmodes of X/T, on which any function on X/T can be developed, 
respect the boundary conditions imposed by the topology. That is, the eigenmodes of X/T correspond precisely to 
those eigenmodes of X that are invariant under the action of the holonomy group F. Thus any linear combination of 
such eigenmodes will satisfy, by construction, the required boundary conditions. In this way we visualize the space of 
eigenmodes of X/T as a subspace of the space of eigenmodes of X, namely the subspace that is invariant under the 
action of F. The computational challenge is to find this F-invariant subspace and construct an orthonormal basis for 
it. In the case of flat manifolds the eigenmodes of X/T — R^/F can be found analytically. In the case of hyperbolic 
manifolds, many numerical investigations of the eigenmodes of X/T = H^/F have been performed 18, 37, 38, 3^ 
113, 13 • In the case of spherical manifolds, the eigenmodes of X/T ~ S^/F have been found analytically for lens and 
prism spaces '36| and otherwise can be found numerically . 

r -v'l 

In the following, we will develop the eigenmodes of X/T on the basis 3^^^^ of the eigenmodes of the universal 
covering space as 

=0 m=-e 



fri s 

so that all the topological information is encoded in the coefficients Cklm' where s labels the various eigenmodes 
sharing the same eigenvalue — k^, both of which are discrete numbers . The sum over i runs from to infinity if 

the universal covering space is non-compact (i.e., hyperbolic or Euclidean). The Cfc^m coefficients can be determined 
analytically for Euclidean manifolds. In the case of a spherical manifold, the modes are discrete, 

k^{iy + l)y/K (5) 



^ We work here in comoving units, but when spatial section are not flat, the curvature of the comoving space is not normalized: it is not 
assumed to be +1 in the spherical case nor is it assumed to be —1 in the hyperbolic case. Thus our eigenvalues may differ numerically 
from those found in the mathematical literature, although of course they agree up to a fixed constant multiple 

* The spectrum is discrete when the space is compact, e.g. the torus or any spherical space. In a non-compact multi-connected space 
such as a cylinder it will have a continuous component. 
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where v is a nonnegative integer (the cases i/ = and u = \ correspond to pure gauge modes [SjI), so that k\ = 
viy + 2)K, and the sum over I is finite and runs from to since the multiphcity of a mode k is at most its 
multiphcity in the universal cover, which is {v + 1)^. Among spherical spaces, the case of prism and lens spaces 
are the simplest since one can determine these coefficients analytically . The worst situation is that of compact 
hyperbolic manifolds, which is analogous to the Euclidean case since the universal covering R"^ is not compact but 
for which no a naly tical forms are known for the eigenmodes. One then needs to rely on numerical computations (see, 
e.g., Refs. IHIiJ)- 

Our preceding works provided a full classification of spherical manifolds (3 St] and developed methods to compute 
the eigenmodes of the Laplacian in them Among all spherical manifolds, we were able to obtain analytically the 
spectrum of the Laplacian for lens and prism spaces |36l | which form two countable families of spherical manifolds. 
The goal of the present article is to simulate CMB maps for these two families of spherical topologies, as well as for 
the Euclidean topologies. 

We first review briefly the physics of CMB anisotropies (Sec. mainly to explain how to take into account the 
topology fSec. Ill Bll once the coefficients ^^ll^ are determined. We then detail the computation of these coefficients, 
focusing on the cases where it can be performed analytically, that is for flat spaces and lens and prism spaces. We 
then present results of numerical simulations fSec. lIV|l as well as simulated maps. We discuss these maps qualitatively 
and confirm that the expected topological correlations (namely matching circles Q) are indeed present. The effects 
of the integrated Sachs- Wolfe and Doppler terms, as well as the thickness of the last scattering surface, are discussed 
in order to give some insight into the possible detectability of these correlations. Figure ^ summarizes the different 
independent steps of the computation as well as their interplay. 

Notations 

The local geometry of the universe is described by a Friedmann-Lemaitre metric 

ds2 ^ -c^di' + {t) [dx' + s\ (x)dr!'] , (6) 

where a{t) is the scale factor, t the cosmic time, = dO^ +sin^ Odip^ the infinitesimal solid angle, x is the comoving 
radial distance, and where 

{sm\i(yJ\K\x)l (hyperbolic case) 

X (flat case) (7) 

sin(\/^ x) /VK (spherical case) 

In the case of spherical and hyperbolic spatial sections, we also introduce the dimcnsionlcss coordinate 

X = VWW. (8) 
which expresses radial distances in units of the curvature radius. 



II. CMB ANISOTROPIES 



The equations dictating the evolution of cosmological perturbations are local differential equations, and remain 
unchanged by the topology of the spatial sections. Indeed the goal of this section is not to derive the whole set of 
equations that has to be solved (see, e.g., Refs. ^3)E3)0)E3 for reviews) but to explain the steps that must be 
changed to take into account the topology. 



A. Local physics of cosmological perturbations 

The CMB is observed to a high accuracy as black-body radiation with temperature To = 2.726 ± 0.002 K [46l| . 
almost independently of the direction. After accounting for the peculiar motion of the Sun and Earth, the CMB has 
tiny temperature fluctuations of order ST/Tq ~ 10~^ that are usually decomposed in terms of spherical harmonics 

s-rp CO i 

jr{d,v) = Y. E «?™^r(e,¥'). (9) 
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I: Universal cover (X) eigenmodes 
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II: CMB in X 
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(see Sec. lII Afl 
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V: Output: aim, maps, C'e. 



(see Sec. Eg 



III: Eigenmodes of X/T 
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IV: CMB in X/T 
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[X] (r^[T] 
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Perturbation equations 
(local physics, unchanged) 



Initial conditions 
(unchanged) 

(see SecCHJ 



FIG. 1: The computation of CMB anisotropics in multi-connected spaces follows a series of independent steps. From the 
knowledge of the spatial curvature K, one determines the universal covering space, X, as well as the eigenmodes of the 
Laplacian in this simply-connected space (upper left box); these functions are well-known and frequently used in standard 
cosmology computations, and are recalled in Appendix^ [we have introduced lo — k/ yJ\K\ for the hyperbolic case and 
V = k/^/K — 1 for the spherical case [recall Eq. Q], and used x defined by Eq. Q]. Once some cosmological parameters and 
a scenario of structure formation has been chosen (upper right box), the CMB anisotropies in the universal covering space can 
be computed. This step is also a standard step that can be achieved numerically by a number of codes (see Sec. Ill At . An 
independent computation (lower left box) is the determination of the eigenmodes of the Laplacian compatible with the topology 
of space, specified by the choice of the holonomy group F. Our approach is to encode all the topological information in a set 
of parameters CfcJm- Their computation is described in Sec. 1111 B 21 and can be performed either numerically or analytically 
according to the case at hand. The implementation of the topology in a standard CMB code (lower right box) is described in 
Sec. Ill Bl and yields the complete correlation matrix of the atm from which one can (center right box) compute Ci, simulate 
maps, etc. 
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This relation can be inverted by using the orthonormahty of the spherical harmonics to get 

ai^^ J ^Y;''*smededip. (10) 

The coefficients agm obviously satisfy the conjugation relation 

= (-ira,*„. (11) 

The angular correlation function of these temperature anisotropics is observed on a 2-sphere around us and can be 
decomposed on a basis of the Legendre polynomials Pi as 

f^(7i)f^(72)) = C°^'i9i2) ^^J2(^i+l)Cf^Pdco8e,,), (12) 

' 7i.72=cosei2 I 

where the brackets stand for an average on the sky, i.e., on all pairs of directions (71,72) subtending an angle 9i2- 
The coefficients C^^^ of the development of C°^^{di2) on the Legendre polynomials are thus given by 

These C^^^ can be seen as estimators of the variance of the a^™^ and represent the rotationally invariant angular 
power spectrum. They have therefore to be compared to the values Ce predicted by a given cosmological model, which 
is specified by (i) a model of structure formation which fixes the initial conditions for the perturbation (e.g., inflation, 
topological defects, etc), (ii) the geometry and matter content of the universe (via the cosmological parameters) and 
(in) the topology of the universe. 

For the case of a simply connected topology, the Ci are usually computed as follows. First, one starts from a set 
of initial conditions given by an early universe scenario. Typically these are the initial conditions predicted in the 
framework of an inflationary model, although this is of no importance in the discussion that follows. Second, the 
modes of cosmological interest are evolved from an epoch before they enter into the Hubble radius till now. Third, 
one computes the variance for the ae m ■ The details of this procedure are well known [4^ lisl |3 . 

For simplicity let us sketch the case of a flat universe. The temperature fluctuation in a given direction of the sky 
can be related to (i) the eigenmodes exp(ik • x) of the Laplacian (here x represents a vector in the usual Cartesian 

coordinate system) by a linear convolution operator Oj. (exp(ik • x)) depending on the modulus k only [as well as 
on the universal cover (here, R^) and the cosmological parameters], and (ii) a three-dimensional variable Ck related 
to the initial conditions, by the formula 

where 7^0 (fc) is the gravitational initial power spectrum, normalized so that 7^0 (fc) oc k^^ for a Harrison-Zel'dovich 
spectrum, and where in most inflationary models the random variable et describes a Gaussian random field and 
satisfies 

{euel,)^S^{k-k') (15) 

with e_k — when the space is simply connected. This relation stems from the fact that the temperature fluctuation 
is real, but it may be expressed differently for multi-connected spaces. Decomposing the exponential by means of 
Eq. I|C7|I and using Eq. I|C5I) allows to rewrite the temperature fluctuation as 



6T 
T ^ 



lip) = E^'/^'dfe^^^of M y|i,(fcx)j J dn^Yr{Ok,^^)ek 



j fc^dfe^^of 1 (yf^l) eUk), (16) 

i.m 



^ One can show that this is indeed the best estimator of their variance when the fluctuations are isotropic and Gaussian l47f . 



7 



where we have defined 

e,™(fc) = y"dfikrr*(^k,<^k)ek, (17) 

which is the "average" of the random fields Ck over all the k of the same modulus. This quantity can therefore be 
identified as a two-dimensional Gaussian random variable satisfying (e^m(fc)e|/„/ (fc')) = S^{k — k')5tg'5mm' /k'^ ■ It 
follows that the coefRcients aim take the general form 



rdfc 



^V^{k)Gi{k)eira{k), (18) 



with 



G,{k) = of '1 (4f 1) , (19) 
and here the function Gf{k) can be approximated by (see, e.g., [islls^ ) 

Geik) = MHvo-V.ss)) (^^(fc,ry,ss) + <i>(fc,%ss) + *(fc,%ss)) + - %ss))^^^^^^^ 

J^(fc(?7o-?/))($(fc,?7) + *(fc,r;)) d?7, (20) 

which is indeed a linear convolution operator acting on i?).^ as announced above. Here, t^lss ^^nd 770 are the conformal 
time at the last scattering epoch and today, respectively, ji is the spherical Bessel function of index £, $ and 5* are 
the two Bardeen potentials, and Vh is the velocity divergence of the baryons. The only modification of note when one 
considers a non flat universe is that the ji are to be replaced by their analog for non flat geometries, the so-called 
ultraspherical Bessel functions , and for a closed universe the integral over k is replaced by a discrete sum. 

In conclusion, and without loss of generality, the temperature fluctuation can be decomposed as in Eq. (|16() whatever 
the curvature of space. For a simply connected topology and Gaussian initial conditions the addition property l|('5|) 
of the spherical harmonics imposes that 

{aima}, m') = C^i^ei'^rnm' , (21) 
and therefore the Cg coefficients encode all the information regarding the CMB anisotropics. 



B. Implementing the topology 

The topology does not affect local physics, so the equations describing the evolution of the cosmological perturba- 
tions are left unchanged. As a consequence, quantities such as the Bardeen potentials $, ^, etc., are computed in the 
same way as described above, and the operator 0[, is therefore the same. However, a change of topology translates 

into a change of the modes that can exist in the universe. In particular, the functions 3^j,£„j are typically not well 
defined on the quotient space X/T. Therefore, the only change that has to be performed is the substitution 

yl^L^^fl (22) 

rpi 

where the j form an orthonormal basis for the space of eigenmodes of the Laplacian on the given topology X/T. One 
must then remember that the mode functions Tj. j can be decomposed uniquely by Eq. Q and that the convolution 
operator o["^' is linear. When the multi-connected topology is compact, it follows that Eq. (|16|l will take the form 

where now Ck is a three-dimensional random variable which is related to the discrete mode k. Equivalently one can 
write ik = iks where k is the modulus of k and the index s describes all the eigenmodes of the Laplacian for fixed 
modulus k in the topology X/T of volume V. These random variables satisfy the normalization 



(24) 
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For a given value of the wavenumber k, there are fewer eigenmodes in the multi-connected case, so that s has to be 
seen as a "subset" of the set {i,m}. 

By inserting the expansion of T^'"] in terms of the covering space eigenmodes, as given by Eq. Q), we obtain 

f (^'^) = ^EE^O^oi"' {yi"l) v^^.- (25) 

k,s i.ni 

It follows that the a^m, seen as random variables, are given by 

fc s 

Note that the sum over s is analogous to the sum over angles defining the two-dimensional random variable eim in 
Eq. (|16|l . Since the aim are linear functions of the initial three-dimensional random variables, they are still Gaussian 
distributed but they are not independent anymore (as explained before, this is the consequence of the breakdown of 
global isotropy and/or homogeneity). The correlation between the coefficients agm is given by 

fc s 

Clearly these correlations can have non-zero off-diagonal terms, reflecting the global anisotropy induced by the multi- 
connected topology, so that Eq. (|21|l no longer holds and the observational consequences of a given topology on the 
CMB anisotropics are given by the correlation matrix QJ. This means in particular that for fixed £, the a^™ might 
not have the same variance, although they all follow Gaussian statistics as long as the initial conditions do. This 
translates into an apparent non Gaussianity in the sense that the Cg will not follow the usual distribution. Strictly 
speaking, this is not a signature of non Gaussianity but of anisotropy. 

Note also that the correlation matrix l|27|) is not rotation invariant. It will explicitly depend on the orientation of 
the manifold with respect of the coordinate system. However, knowing how the spherical harmonics transform under 
a rotation allows us to compute the correlation matrix under any other orientation of the coordinate system. To finish 
let us note that one can define the usual Ce coefficients in any topology by the formula 

= ^g-^ E '^^^^ ' (28) 

m 

which is easily shown to be rotationally invariant. 



III. EIGENMODES OF MULTI-CONNECTED SPACES 



A. Flat spaces 

fri s 

The purpose of this subsection is to compute in detail the coefficients Ckim '-'^ ^ cubic 3-torus of comoving 

size L, referred to as Ti. The method generalized easily to any compact fiat manifold. 
In the case at hand, the topology implies the "quantization" of the allowed wave vectors 

k = — (nix + n2y + n^z) = — n = fcn, (29) 



with 



(ni,n2,n3), (30) 

Vim, (31) 

{ni,n2,n3)/n, (32) 
2tt 



k EE Vk^^ — n, (33) 
L 

so that the label s can be chosen to be the "unit" integer vector fi. The label multiplicity mult(fc) = cardjs} of a 
mode k is given in this case by the number of representations of by 3 squares, i.e., by card{n} (see Fig. |2below). 
The corresponding normalized eigenmodes in Cartesian coordinates are thus simply given by 
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with k given by Eq. H29|l . Using the decomposition (|C7|) of the exponential and plugging in the closure relation ljC5p . 
one gets that 



Skim — ?27™/i^m " ' V"^ 

where n can also be defined by the two spherical angles (0n, fn) which are explicitly given by 



tan^n 



tan (ySn = 



n2 
ni ' 



(35) 

(36) 
(37) 



This expression could also have been obtained by simply considering the Fourier transform of y\^i^ as given by 
Eq. (jC6p . One can check that the normahzation of the basis T^f^' [i.e., J Tj.^'T{J^|,*d^x = J°(k— k')] implies that 



the coefficients Cfc^^IrT satisfy the closure relation 

oo t 

E E d'^i'^d^y™ * = -J^lcOS^n - COS0„O5°(¥'n " ^n')4fc'. 

1=0 m=-t 

The Dirac distribution in the above expression can be shown, by using either Eqs. (jC4IC5|) or Eq. 

21+1 



5°(COS0„ - COS0„O<5°(¥'n - <P„') = E ^^^^(fi-fi' 



£=0 



^From these results, we deduce that the correlation matrix of the aim is given by 



f^) f ^) [ ^) Err*(n)r/(n) 



n \ / \ / \ / 

Using Eq. (jC5|) and the fact that ^ — the Ci coefficients are simply 



(38) 

alone, to be 
(39) 

(40) 



IT 



2Tm 



(41) 



which was used in many earlier works [^^IHI^ on the influence of topology on the CMB. 
Note that spherical harmonics satisfy the following symmetry relation 

Yp{e, -ip) = e-2""^r^™(0, if). 



(42) 



Due to the symmetry of the torus with respect to the y — Q plane, in the sum over n in Eq. H4U|) a term (^nj'/'n) 
will always be associated to a term ((?„, —Lp„) leading to a term of the form 4 cos {{m! — m)Lpn) Yl"'*{9n, 0)Yf™ (^n, 0), 
the only exception being the term arising when n2 = 0, which is real. ^From this result, one easily shows that the 
correlation matrix satisfies 



Along similar lines, the relations 



Yf'"(6',7r) = c*™'^V(6',0), 



imply 



so that 



C 



l + (-l)'"-™ ! + (-!) 



<e' m' 

£ m ' 



(43) 



(44) 
(45) 



(46) 



Cf„™V0 ^ m-m' = 0mod(2) and £ - £' = mod(2). 



(47) 
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Furthermore Eas. H43IC2|l imply 



cfr' = c^fjr'- (48) 



Let us emphasize that these properties of the correlation matrix still hold even if the torus is not cubic. However, a 
cubic torus is invariant under a 7r/2-rotation about the z axis, so if (rii,n2,"-3) corresponds to a wavenumber then so 
does (n2, — n-i, Ha), and one has 

CfrVo ^ m-m' = mod(4). (49) 



B. Spherical spaces 

The goal of this section is to recall the basic analytical results concerning the lens and prism spaces fSec. 1111 fj 2ll . 
In these spaces, the eigenmodes and eigenvalues of the Laplacian operator can be determined analytically using 
toroidal coordinates [361] . CMB computations use spherical coordinates, so we must perform a change of coordinates 
and a change of basis (detailed in Appendix 0. Fortunately, this can also be achieved analytically to compute the 
coefficients 



1. Generalities 

In our preceding article |3.5| . we presented in a pedestrian way the complete classification of three-dimensional 
spherical topologies and we described how to compute their holonomy transformations. 

The isometry group of the 3-sphere is SO (4). Every isometry in SO (4) can be decomposed as the product of a 
right-handed and a left-handed Clifford translation, and the factorization is unique up to simultaneous multiplication 
of both factors by -1. Furthermore, the space itself enjoys a group structure as the set of unit length quaternions. 
Each right-handed (resp. left-handed) Clifford translation corresponds to left (resp. right) quaternion multiplication 
of , so the group of right-handed (resp. left-handed) Clifford translations is isomorphic to . It follows that S0(4) 
is isomorphic to x 5^/{±(l, 1)} and thus the classification of the subgroups of S0(4) can be deduced from the 
classification of subgroups of . There is a two-to-one homomorphism from to SO (3); the finite subgroups of 
SO (3) are the cyclic, dihedral, tetrahedral, octahedral and icosahedral groups, so the finite subgroups of are their 
lifts, namely 

• the cyclic groups Z„ of order n, 

• the binary dihedral groups Z?*„ of order 4to, m > 2, 

• the binary tetrahedral group T* of order 24, 

• the binary octahedral group O* of order 48, 

• the binary icosahedral group /* of order 120, 

where a binary group is the two-fold cover of the corresponding plain group. 

^From this classification, it can be shown that there are three categories of spherical 3-manifolds. 

• The single action manifolds are those for which a subgroup R (resp. L) of acts as pure right-handed (resp. 
pure left-handed) Clifford translations. They are thus the simplest spherical manifolds and can all be written 

as s3/r with r = Zn,D;^,T*,o*,r. 

• The double-action manifolds are those for which subgroups R and L of act simultaneously as right- and 
left-handed Clifford translations, and every element of R occurs with every element of L. There are obtained for 
the groups F = Fi X Fa with (Fi,F2) = (Z™,Z„), {D*^,Zn), (T*,Z„), (0*,Z„), (r,Z„), with gcd(m,n) = 1, 
gcd(4m,n) = 1, gcd(24,n) — 1, gcd(48,n) = 1, and gcd(120, n) = 1, respectively. 

• The linked-action manifolds are similar to the double action manifolds, except that each element of R occurs 
with only some of the elements of L. 



The classification of these manifolds is summarized in Fig. 8 of Ref. |3^ 
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2. Lens and prism spaces 

In this article, we focus on prism spaces S^/D^^ and lens spaces L{p, q) . The latter are defined by identifying the 
lower surface of a lens-shaped solid to the upper surface with a rotation of angle 2-Kq/p for p and q relatively prime 
integers with Q < q < p. Furthermore, we may restrict our attention to < g < p/2 because for values of q in the 
range p/2 < q < p the twist 2'iTq/p is the same as — 27r(p — q)/p, thus L(p,g) is the mirror image of L{p,p— q). Lens 
spaces can be single action, double action, or linked action; Fig. 9 of Ref. [33 summarizes their classification. 

The eigenmodes and eigenvalues of prism and lens spaces can be obtained analytically by working in toroidal 
coordinates '36'|. Starting from Cartesian coordinates, in which the equation for the 3-sphere is + + + w'^ = 1, 
the toroidal coordinates {xtjOt^^Pt) are defined via the equations 

X = cosxt cos6't, (50) 

y = cosxt sin^T, (51) 

z — sin Xt cos (y^T , (52) 

w — sinxx sin(^T, (53) 

with 

< Xt < 7:/2, (54) 

< < 2n, (55) 

< ip-r< 27T. (56) 

Ref. gives the eigenmodes of S'^ explicitly as 

Q.^T^T =5.^TmTC0sl^-lxT sinl™-lxTPi'""'^""'Hcos2xT)/(KT|eT)/(|mT|</5T), (57) 

where p is the integer parameterizing k = {v + l)v^ as in Eqn. p^'™"^' ' l^'^'-' is the Jacobi polynomial, and / 
stands for the cosine (resp. sine) function when i-^. or ttit > (resp. £t or ttit < 0). For each value of the indices 
and ttzt range over all integers satisfying 

I^tI + |"It| < I' and (58) 
\e^\ + |tot| = mod(2), (59) 



and for convenience we define 



d=^ii^-\U\-\m^\). (60) 



The normalization coefficients £^ m-r are given by' 
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_ g^T^^T l 2{j^ + l)dl{\e.r\ + \m-r\ + d)l 

^{v + l)\j (|£,|+d)!(K|+d)! ^""'^ 

with Gi = 1/ y/2 if i = and (Xi = 1 otherwise. 

Using these definitions, Ref. [sBl shows that for lens spaces the explicit set of coefficients 'nfj^rm ^^'^ that 

^fel = X! ^rjTmT2^^TmT(XT,6'T,</3T), (62) 



can be obtained as follows: 



Note the factor + 1) which differs from Ref. l3d due to a different choice of normahzation. 
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Theorem 1: lens spaces 



The eigenspace of the Laplacian on the lens space L{p,q) has an orthonormal basis that, when lifted to Zp-invariant 
eigenmodes of the 3-sphere, comprises those eigenmodes in the left column for which the corresponding condition in 
the right column is satisfied, subject to the restriction that an eigenmode rm exists if and only if the integers v , 
£tj Ound rriT satisfy |-^t| + {"titI < and \£t-\ + IttitI = ly mod(2). 



basis vectors 








condition 


Qi/OO 








always 


QuIt 








= mod(p) 


Qi/ niT 








qniT- = mod(p) 


[Qvl-T m-T + Qv -£t 




( Qi/ —^T■ *TiT 


Qul-T -mx)/ ^2 


It = qm^ mod(p) 


(Qvl-T m-T ~ Qv -It 


-mT)/\/2i 




Qvi-T -mx)/ V2 


It = —qm^ mod(p) 



An analogous theorem was demonstrated for prism spaces and can be found in Ref. 36]. 

Unfortunately, for practical purposes the eigenmodes of the lens and prism spaces are needed in spherical coordi- 
nates, while they are most easily obtained in toroidal coordinates. As explained in the Introduction, one needs the 

coefficients ^if^^j of the decomposition Since y^^im ''^^'^ QviTim are two orthogonal bases of dimension {v + 1)^, 
all of whose elements have the same norm, there is an orthogonal transformation taking one to the other 

£,m 

The "transpose" of this transformation a takes a given eigenmode's Q-based coefficients 77 to its 3^-based coefficients ^: 

(64) 



E\ ^ [r] u 

/ . Civ ImlTm-T^stT-n 



It = — u TriT^ — v 

The orthonormality of the basis tJJ"] implies that the coefficients satisfy 



(65) 



£=0 m=-e 



This relation is simpler than the closure relation H38|) obtained in the flat case because, for a given fc, the space of 
modes is finite-dimensional. The computation of the coefhcients ai, i ,„ £t mx appears in Appendix ^ Because both 
the Q basis and the y basis are orthonormal, the transformation a is orthogonal: 



(66) 



where ev{lT,mT) = 1 if the conditions H58I59|I are satisfied and otherwise. 

With these coefficients, the CMB computation goes as in the flat case, except for the fact that some integrals have 
to be replaced by discrete sums. One easily gets that 



so that 



and 



(27r 



,3 °° 



V 



v=2 s 



(2^ 
V 



3 °o 



(67) 



(68) 



(69) 
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as first obtained in Ref. Note also that, following Refs. a scale invariant spectrum will in that case be 

defined as 

To finish, let us discuss the properties of the random variable ek- Since the eigenmodes in toroidal coordinates, 

real valued, it follows from Eq. (C2) that 

'^vlmlTm.T ~ {^'^)™Ciyl-mlTmT- (71) 



QuiT-rriTJ ^iiti the coefficients fyU^^^^ are real valued, it follows from Eq. (C2) that 



It follows that 



f[r]s * _ c_i v'y'^^'* (72) 



[ri 

whatever s and thus that the eigenmodes TJ. ^ are real- valued. This implies that Ck is a real random variable, contrary 
to the preceding example of the torus. 



IV. NUMERICAL COMPUTATIONS 



A. Implementation 

The correlation matrix for £ < i'max has ^^ax coefficients. However, the parity and symmetry relations 

/-it' m' / -1 \m — m'y^i'' — m' * ^'7Q^ 

^tm — \^^) '^e.-m I'-JJ 

r^i' ra' r~ftra * (nA\ 

reduce the problem to computing only a quarter of them. Then, for a given topology, symmetries can further reduce 
the number of coefficients to compute. For example, with a cubic torus, Eqs. H47I49(I insure that only one coefficient 
out of eight is nonzero, and the symmetries (|48|l also give the coefficients when one changes the sign of both m and 
m! . This leaves only £j^ax/64 coefficients to compute. For example, a COBE scale map (^max ^ 30) requires 12500 
coefficients, while a Planck scale map (£max ~ 1500) requires ^ 8 x 10^° coefficients. 

Each coefficients is computed using Eq. H40|) , which involves a sum over all the wavemodes k. For a given resolution 
^max the modulus of the largest wavemode is given by fcmax 3£niax/?yo- Moreover, the density of wavemodes is 
proportional to the size of the torus, so that we have 0(l\^^^L^) modes. Therefore, the computational time and the 
memory requirement scale as P^^^L? and (^^^L^ ^ respectively. This is obviously a serious limitation of our algorithm. 
For example, computing the correlation matrix for a COBE scale map on a relatively small torus (L — 2Rh, where 
Rh is the Hubble radius) takes around 10 hours on a 900 MHz CPU and allocates 60 Megabytes of memory. 

When the topology is simply connected, it is a well-known fact that the Ci are in general a smooth function of the 
multipole £. This reduces the computational time because for ^max = 1500 one needs to compute only ~ 50 coefficients. 
For the topologies we have studied, we did not find any evidence for a smooth structure of the correlation matrix, at 
least at the relatively large scales we considered. At which scale one can reliably approximate the correlation matrix 
by its isotropic diagonal part (the Ci) remains an open question. 

Also, if one wants to simulate CMB maps from the correlation matrix, one needs to diagonalize it. This procedure 
can also take a lot of time because it is an ^^^x process. For the case of the torus, however, this problem is not serious 
as the symmetries of the torus insure that the matrix is block diagonal, with eight blocks if the torus is cubic or four 
blocks otherwise. 

Strictly speaking, one does not need the correlation matrix to compute maps. One can do it directly by using 
Eq. H18|) . This amounts to performing a realization of the three-dimensional random field describing the cosmological 
perturbations, and projecting it onto the sphere. In this case, one has only ^max coefficients to compute (the aim) 
instead of the correlation matrix, so that the memory requirements are roughly the same (one only needs to store the 
value of the random field for each mode), but the computational time scales as i^^^L^. In this case, computing maps 
for a cubic torus of size L = 2Rh till Anax = 120 takes 3 hours on a 1.7 GHz CPU and allocates 300 Megabytes of 
memory. 

[ri s 

In the case of spherical spaces, the coefficients Ckim must also be computed numerically. This involves determin- 
ing both the coefficients rj^si^mT ^^'^ af^m^TwiT- This computation can be reduced by taking into account their 
symmetries, as described in Appendix IbI which imply that f — |m| -I- \£t\ is odd, \m\ = |tot|, Kt| + |wt| < t^, and 
\£t\ + I^tI = mod(2). The computation can be performed analytically with the use of symbolic computation 
software such as Mathematica. In the case of the lens space L{17, 5), the computation up to = 43 and = 55 takes 
3 and 12 hours respectively on a 1 Ghz CPU with a negligible amount of memory. 
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FIG. 2: Multiplicity of the mode fc as a function of n (cubic torus, top panels) or f (lens space L{p,q), bottom panels). For 
the cubic torus, n is of the form \/]V, where A'^ is an integer, and for lens spaces (and more generally spherical spaces), u is an 
integer. Left panels show the multiplicity of each mode for a given value of either n of v. We have also given an estimate of the 
"average" number of modes, given by 27rn for the torus and {u+ Vf /p for the lens spaces, respectively. For lens space L(17, 5), 
the mode multiplicity closely follows the analytical estimate, whereas for lens space 1/(12, 1), modes exist only for even values 
of u. For the cubic torus, the mode multiplicity is even more irregular, and varies between and ~ 20n. For example, it is 
always when n — V8m + 7. Right panels show the integrated mode multiplicity, that is the number of modes smaller or equal 
to some value (2-Kn/L and {y + 1)^/K, respectively). Here, the anal ytic al estimates {Airn? /Z and {v + !)(;/ + 2)(2i^ + 3)/6p, 
respectively) provide a much better estimation (see Section 5 of Ref. [lj| for a detailed discussion of these properties). 



B. Expected results 

The three main effects that are expected on a CMB map computed in a multi-connected topology are (i) the 
appearance of l-t' and m-m' correlations reflecting the breakdown of global isotropy, (ii) the existence of a cutoff 
in the CMB angular power spectrum on large angular scales (low £) , and (iii) the existence of pattern correlations 
such as pairs of circles where the temperature fluctuations are strongly correlated as they represent the intersection of 
the last scattering surface with itself and therefore show the temperature of the same emission region from different 
directions. Note that the effects (ii) and (iii) will show up only if the topological scale is smaller than the radius of the 
last scattering surface while the effect (i) may be present even if the topological scale is a bit larger than the diameter 
of the last scattering surface. 

So far, the main constraints that have been given on multi-connected topologies come from the absence of a cutoff 
at large angular scales in the COBE spectrum. This gives strong constraints on the minimal size of the topology as 
the cutoff is given by the angular size of the torus projected on the last scattering surface. However, as previously 
discussed, this cutoff in the "true" temperature fluctuations can be compensated, at least partially, by an integrated 
Sachs- Wolfe effect which arises, for example, when the cosmological constant is large. 

The third, and up-to-now never computed, effect of a multi-connected topology is the appearance of pairs of circles 
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which are correlated in temperature. This correlation, however, is not perfect. It would be perfect if the temperature 
fluctuation were a pure scalar function on the last scattering surface 2-sphere around us which would be the case only 
if (i) the temperature anisotropics were given only by the Sachs- Wolfe effect [first term of Eq. (jSOJ] and (ii) the last 
scattering surface were infinitely thin. 

It is well-known that the temperature fluctuations observed in a given direction are in fact a combination of several 
effects: first, one has the intrinsic temperature fluctuations of the emitting region, which is eventually affected by a 
gravitational redshift. These two contributions form the so-called Sachs- Wolfe effect [first term in the righthand side 
of Eq. H2U|) ]. Second, if the emission region is not at rest with respect to the observer, one will observe some apparent 
temperature fluctuations which in fact result from a Doppler shift [second term in the righthand side of Eq. (|20() ]. 
Third, several events can alter the photons energy and trajectory while travelling toward us. In particular, they can 
be slightly disturbed from their trajectory (lensing) and, more importantly, they can exchange energy when they cross 
time- varying potential wells. This last effect is usually referred to as the integrated Sachs- Wolfe effect [third term in 
the righthand side of Eq. see Fig. |2| below]. Obviously, the Sachs- Wolfe effect is a scalar quantity that depends 
only on the emission region. Therefore, it should be the same whatever the direction of observation. By contrast, the 
Doppler effect will explicitly depend on the direction of observation. If one observes two directions which correspond 
to the same point of the last scattering surface and which form a small angle, then one expects that the Doppler 
contribution will be almost the same. If the matching points are 90 degrees from each other, then one expects on 
average no correlation at all, whereas the Doppler effect between two antipodal points will become anti-correlated. 
Finally, since photons originating from the same emission region but observed from different directions will travel 
through different regions of space, they will undergo different integrated Sachs- Wolfe effects, so that no significant 
correlations are expected from this effect, which is therefore considered a noise term for our purposes. 

Actually, one aim of this work is precisely to compute the typical amount of correlation one can expect on pairs 
of circles. Note that this correlation is likely to depend on scale: on large scales, one should be annoyed by the late 
integrated Sachs- Wolfe effect; between the Sachs- Wolfe plateau and the first Doppler peak (and, at a lesser extent, 
at every dip between two Doppler peaks), the Doppler effect dominates; at the first peak, there is usually (especially 
when the matter content is low) a significant contribution of the early integrated Sachs- Wolfe effect (see Fig. I^l ; at 
very small scales, one feels the finite width of the last scattering surface (see below), etc. Also, as explained above, the 
relative position of the circles will play a role because of the Doppler contribution. It is therefore interesting to look 
at the best way to find matching circles on a realistic CMB map. We leave this important point to future work 51]. 

The second (and probably less important) effect that reduces the correlation between the circles is the finite width 
of the last scattering surface. As far as we know, this effect has not yet been carefully analyzed. It plays a role when 
one looks at fluctuations on scales smaller than the projected width of the last scattering surface. In this case when 
looking in a given direction, one picks up fluctuations which are situated "on one side" of the last scattering surface, 
but for pairs of circles, one sees opposite sides of the last scattering surface. On larger scales, the effect is negligible 
as one averages temperature fluctuations on regions much larger than the thickness of the last scattering surface. 

V. RESULTS 

We now outline some of the results we have already obtained from our simulations. The main aim of this Section 
is to provide a series of tests to check our simulations. A more detailed analysis of the structure of the correlation 
matrix C^^" as well as a search for accurate tests to detect the topology are left for future work [sif . 

A. Flat case: cubic torus 

In all the simulations we performed, we have considered a flat ACDM model with Ha = 0.7, a Hubble parameter of 
Hq = 100/ikms~^ Mpc~^ with h = 0.62, a baryon density Wb = ^hh^ = 0.019 and a spectral index ng = 1. With this 
choice of cosmological parameters, the Hubble radius is Rh ~ 4.8 Gpc, the "horizon" radius (under the hypothesis 
of a radiation dominated universe at early times) is Rh ~ 15.6 Gpc, and the radius of the last scattering surface is 
i?Lss ^ 15.3 Gpc. The volume of the observable universe is therefore Vobs ~ 15 x 10'^ Gpc'^. 

Let us first compare the Cg in the simply connected topology to the Ce in a torus (Fig. As expected, we see 
a cutoff at some angular size which corresponds to the angular size of the torus on the last scattering surface. This 
corresponds to the multipole 



(75) 



16 




10 100 
I 

FIG. 3: CMB anisotropies in the simply connected (i.e., usual) topology and in toroidal universes based on cubic fundamental 
domains of various sizes, expressed here in units of the Hubble radius. With our choice of cosmological parameters, the situation 
L = 8 corresponds to a torus larger than the observable universe, which shows a small depletion of power on large scales. For 
smaller tori, the cutoff is much sharper. 



where L is the length of the cubic torus' fundamental domain |49l |50j. Note that even when the torus is larger 
than the size of the observable universe, the spectrum exhibits a loss of power on large scales. This is because the 
Harrison-Zel'dovich spectrum exhibits a significant amount of power at large scales (by definition, it is scale invariant), 
and in practice, the modes that contribute to the quadrupole of the CMB anisotropies can be as large as ten times 
the size of the observable universe (the exact number depends mostly on the spectral index and on the amplitude of 
the integrated Sachs- Wolfe term). Therefore, this leaves hope to detect the topology "beyond the horizon" where the 
circles method would fail. 

It is not easy to predict the amplitude of the power at scales larger than the cutoff because it depends mostly on the 
amplitude of the integrated Sachs- Wolfe effect, which is difficult to estimate even when the topology is multiconnected. 
Another consequence of a multiconnected topology are oscillations in the spectrum. These come both from the fact 
that there is a sharp cutoff in the spectrum (which causes oscillations in Fourier/Legendre space) and that the 
spectrum is "spiky" on large scales. Should we consider a simply connected universe with a cutoff at some scales, 
then the corresponding Ci would be less irregular. Finally, note that on small angular scales, the spectrum tends to 
behave as in the simply connected case, but computing this scale remains an open problem at the moment. 

So far, we have considered only the Cg, which represent only some average of the diagonal part of the correlation 
matrix. The true diagonal part of the correlation matrix is given by the C^™ which represent the variance of the 
aim- An example of their behavior is shown in Fig. 0] for the same topologies as in Fig. |31 Several features appear on 
this figure. First and most importantly, the dispersion in the variance of the a^m at fixed t is very large. It appears 
that it is maximal at the cutoff scale -^c, where the dispersion in the variance of the can be as large as two 
orders of magnitude. This dispersion slowly decays at larger multipoles, where one "tends" (in the sense of observable 
quantities) towards the simply connected case, and surprisingly also decays at scales larger than the cutoff. With the 
hypothesis of a multi-connected universe, this dispersion would be (incorrectly) interpreted as non Gaussianity. Since 
at present no non Gaussianity or anisotropy was observed in the data, this allows new constraints of the size of the 
fundamental domain. 

We did not find any convenient way to represent the off-diagonal terms of the correlation matrix. We therefore 
switch to showing and analyzing some realizations corresponding to the numerically computed correlation matrix. In 
what follows, we have fixed the size of the torus to L = 2Rh- We therefore have 16 copies of the torus in the 
observable universe, and from Eg. 1)75(1 this corresponds to Ic ^ 9, as can be checked in Fig.|21 Although such a model 
is now excluded by the data we analyze it in detail mostly for pedagogical purposes (and also because 

the computing time scales as L^). In the case where the torus, or more generally the fundamental domain, is smaller 
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FIG. 4: Variance of the a^™ for fixed values of £. The average of these values give the Ce coefficients shown in Fig. |3 We 
insist that this does not correspond to a realization of the random variables describing CMB anisotropics, but to the variance 
of the aern- These coefficients are far from sufficient to build maps of CMB anisotropies as they do not include the correlations 
between different atm- 



than the last scattering surface, one expects to see pairs of circles where the temperature is correlated 7]. Seeing 
these circles at their expected position is therefore the most crucial test of the procedure outlined in Sees. IIIIIIIIIVl 
As already announced, the aim here is not to derive a detailed procedure to detect these circles, but to check our 
algorithm and to explore some of the properties of these matching circles. Here, we have 2i?Lss/i = 3.17 ~ -^10. 
One therefore expects to have 61 pairs of circles and 12 pairs of points having correlated temperature''. In order to 
see the circles, it is convenient to show the last scattering surface as a sphere seen from the outside and to look at its 
intersection with itself after a translation of the form L ■ (ni,n2,n3), as shown in Figs.[SHHl 

In the last Section, we pointed out that the correlation would depend on the amount of Sachs- Wolfe, Doppler 
and integrated Sachs- Wolfe effects. The decomposition of the temperature anisotropies both in the simply connected 
topology and in the toroidal cases are shown in Fig.^ Note that we show only the relative amplitude of these effects 
and not their cross correlation. 

We now turn to the correlation between pairs of circles, as introduced in 0. 

If the topology is not known in advance, the relative position between matching circles can be arbitrary, so that in 
general the search for circles is a six parameter problem: two parameters for the center of the first circle, two more 
for the center of the second circle, one for their common radius, and one for their relative phase (i.e. the twist with 
which they are identified). In the case of a torus, the circles sit directly opposite each other on the sky (eliminating 
two parameters) and there is no twist (eliminating another parameter), so the problem reduces to a three parameters 
search®. We are not going to perform such a study, but rather focus on some features of matching circles in a toroidal 
universe. 



The circles correspond to the intersection of the last scattering surface with translates of the form L ■ (ni, ng, ng), where (ni, 712, na) = 
(0,0,1), (0,0,2), (0,0,3), (0,1,1), (0,1,2), (0,2,2), (1,1,1), (1,1,2), and (1,2,2) plus all permutations and sign changes among each 
triplet (rti, 712, na). The pairs of points correspond to the case where the intersection between the last scattering surface and its translate 
reduces to almost a single point, as is the case for (ni, 712, 713) = (0, 1, 3) and its permutations and sign changes. 
® Another way to see that in a torus it's a three parameter search is to visualize the situation in the universal covering space. Place one 
copy of the last scattering surface with its center at the origin, and imagine a translated copy with its center at some point {x,y,z). 
Each choice of {x,y, z) uniquely determines a circle of intersection (assuming < + + < ^lss)' ^'^'^ conversely each pair of 
circles arises from exactly two points {x, y, 2) and {—x, —y, —z) and no others. Thus the point (x, y, z) serves to parameterize the circle 
search. 
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FIG. 5: Four realizations of CMB maps of the whole temperature anisotropies for a cubic torus. The resolution of each map is 
^max = 120 . The last scattering surface is seen from the outside as well as two of its closest topological images after translation 
of L and —L along one axis of the torus. One can check by eye that the temperature fluctuations are well (but not perfectly) 
correlated along matching circles (i.e., along the intersection between the last scattering surface and its neighbors) located at 
latitudes 6 — ±19°. Note that as expected, there are very few fluctuations on scales larger than the size of the torus (which is 
given here by the distance between the circles). 



A simple estimator for the correlation between pairs of circles that are horizontal with respect to the coordinate 
system is obviously 

r(f) f) ) - ^ [ 26(^1, ^)e(g2,y') , 

where we have set Q = ST/T. If the temperature fluctuations on a pair of circles are completely uncorrelated, then 
on average C = 0. If they are completely correlated, then C = 1, and if they are anti correlated, then C = — 1. 
Each aim can be seen either as the theoretical expectation for a given model, or an observed quantity that can be 
measured from the CMB sky. In the first case, it represents a feature that one can expect from a given model, and in 
the second case it represents an estimator of some features predicted by the topology. Here, we shall concentrate on 
the observed C{9i, 62) that we compute from simulated maps, first to check the validity of our procedure to compute 
the correlation matrix C^^' , and second to convince us that it is possible to see the presence of matching circles 
using simple techniques (although we do not pretend that this method is optimal). 

In principle, two matching circles have the same angular diameter, so that only the case 62 = —9i is relevant, 
but we have chosen to leave 62 as a free parameter to sec to what extent uncorrelated circles might happen to seem 
correlated by chance. 

We first show in Fig. ^| a few examples of the observed function C{9i,92) for a simply connected universe. As 
expected, the correlation is quite large when 9i ~ 92 because the circles are near each other and the real space 
correlation function C{9) (the Legendre transform of the Ce) is not when 6* — > 0. With our normalization of 
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FIG. 6: Same as in Fig. |3 but with the Sachs- Wolfe contribution only. These maps have comparatively less small scale power 
than the previous one as ^ = 120 is close to the first dip in the Sachs- Wolfe spectrum, so that one sees the large scales (up to 
the torus angular size) better. Note that the matching between circles is almost perfect here. 



C (01, 6*2), one has 

lim C(0i,02) = l, (77) 

which will of course remain valid when the topology is multi-connected. When the separation between Q\ and 82 is 
large, one can neglect the correlation between the two circles, and the main contribution to C comes from statistical 
fluctuations: it is always possible that two circles exhibit similar temperature patterns by chance. The variance of 
these statistical fluctuations is probably given by the number of independent pixels on the map and therefore by a 
combination of the scale at which the power spectrum is large and of the resolution of the map (here, £max — 30). In 
any case, the amplitude of the largest statistical fluctuations of C(0i, Q-i) gives an idea of the amplitude of the signal 
needed to detect a multiconnected topology^. For the maps we have generated, the correlation reaches 30% for a few 
pairs of falsely matched circles. 

Figs, im 1121 1 131 and 1 141 show contour plots for several realizations of the correlation matrix in a torus universe. The 
torus aligns naturally with the coordinate system, so one expects correlated pairs of circles at 

= -02^± arcsin f ) , (78) 

for each positive integer n such that the arcsine exists. For our choice of cosmological parameters we have i?Lss = 3.17 



The signal threshold could therefore be reduced by performing the same analysis on a higher resolution map, but because the search 
for circles is in general a six-parameter problem, it might be necessary to search low resolution maps first to find likely candidates, and 
then search higher resolution maps to confirm them. 
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FIG. 7: Same as in Fig.|3 but for the Doppler component only. Left panels show as previously the last scattering surface and 
two of its images from the outside, with large circles, which are therefore well correlated (but not as well as for the Sachs- Wolfe 
or total contributions, however). On the right panel, we show smaller matching circles {9 = ±71°), which are more conveniently 
shown from the inside of the last scattering surface. Here, the anticorrelation between circles is obvious. 



and L = 2 (in units of the Hubble radius), giving 

01 = -6I2 = ±18°,±39°,±71°. (79) 

For these values of 9i and 62 , one expects a perfect correlation for the Sachs- Wolfe contribution: 

Csw = 1. (80) 

This formula holds both when one considers Csw as an ensemble average and when one considers a given realization 
of the density field since in both cases it follows from the fact that one sees the same region from different directions. 
These correlations appear clearly in Fig. ^] which considers only the Sachs- Wolfe contribution. In this case one would 
have even expected perfect correlations for the values of 9i and 62 given in H79|l . This is not what we have, but 
the reason for this is easy to understand: imposing 8(6'i,<p) = Q{—0i,(p) in real space induces in Legendre space 
correlations at arbitrary large multipoles i. Here, for computational reasons, we were forced to truncate the correlation 
matrix at a rather low value of £, so the matching is significant but not perfect. It would presumably increase in 
higher resolution maps. 

If one considers the Doppler contribution to the CMB anisotropics, the situation is somewhat different. As an- 
nounced above, the correlation between two circles depends on their relative angle. More precisely, it is given by 

Coop = ( ^^-\2^? ^ / ' (^-'-) 



(ni • n)2 + (n2 • fi)^ 

where fii and 112 are two constant unit vectors spanning an angle \9i —021, and where the brackets denote an average 
over all the directions of the unit vector fi. After some manipulations, one obtains 

Cdop = tan f ^ - . (82) 
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FIG. 8: Eight pairs of matching circles among the 61 existing pairs for a single realization of the density fields. For clarity, the 
orientation of the last scattering surface is the same in all panels. 
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FIG. 9: Decomposition of CMB anisotropies into the Sachs- Wolfe, Doppler and ISW contributions. When the topology is 
simply connected [left panel], the Sachs- Wolfe contribution is dominant at the peaks and usually on the largest scales. The 
Doppler term is dominant only before the first peak, and has a significant contribution between peaks. The ISW term appears 
mostly at the first peak (this is referred to as the early ISW effect) when the radiation-to-matter transition occurs late (typically 
at low h), and at large scales (late ISW effect) when the universe is not matter dominated at z — 0. For standard values of the 
cosmological parameters, it is not dominant. The right panel shows the initial portion of the spectrum {i < 100) for a model 
with the same cosmological parameters but with a toroidal topology. The Sachs- Wolfe contribution shows a sharp cutoff at 
the expected scale. Note that neither the Doppler nor ISW contributions show a similar cutoff. This is due to two different 
projection effects: for the Doppler term, this comes from the fact that it is always negligible on large scales in k space (it scales 
as k^), and when one goes to £ space, the convolution 1201 always transfers some power from the small scales (where there is 
power in k space) to large scales; for the ISW, the presence of power comes from the fact that it is generating long after the 
last scattering epoch, so that it appears on large scales simply because it describes phenomena that occur near us. However, 
since the Sachs- Wolfe contribution is usually largest on large scales, the final spectrum still exhibits a sharp cutoff. 



Again, for the same reason as for the Sachs- Wolfe contribution, this formula holds both if one considers Cdop as 
an ensemble average or if we consider a given realization of the density field. One recovers as expected that the 
correlation is 1, 0, —1 for \9i — = 0°, 90°, 180°, respectively. For the values of the angle given in Eq. H79|l . one 
obtains C = 0.51, 0.07, —0.49, respectively. These are the results that we obtain qualitatively in Fig.^l where no 
correlation at all is seen for the circles at ±39°, and positive (resp. negative) correlation is seen for the circles at ±71° 
(resp. ±18°). 

Finally, the correlations due to the integrated Sachs- Wolfe effect are shown in Fig. ^1 As expected, no particular 
correlation is seen for the values of 9i and 62 of Eq. H79|l . The contour plots are however quite different from those 
of the Sachs- Wolfe and Doppler contributions. The reason is twofold. First, most of the power lies at the smallest 
multipoles. This translates into the fact that the contours are broad in the sense that they do not vary a lot on small 
intervals of 61 and 62- Second, the fact that most of the power is at large scales implies that a very small number of 
modes contributes to it (since we see a finite region of the universe), so there is a large cosmic variance that makes 
the statistical uncertainty very large (thus serendipitously similar temperature patterns on two unrelated circles are 
easily achieved here). 

Combining all the contributions to the CMB anisotropies allows one to simulate realizations of the exact C{6i, O2) as 
shown in Fig. Since the Sachs- Wolfe contribution is dominant, the spikes are still clearly visible at their expected 
positions, but appear less prominent than in Fig.^] As expected, it seems that the circles at 0i,2 — ±39° are slightly 
less correlated than the other two pairs because their Doppler contribution is not correlated, but this deserves a more 
careful analysis. 

B. Spherical case: lens spaces 

Among the spherical spaces, the procedure presented above can be applied most easily to lens and prism spaces, 
because their eigenmodes are known explicitly. The eigenmodes are known analytically in toroidal coordinates (see 
Section llll B|) . and Appendix IbI shows how to convert them to spherical coordinates. In this section, we present some 
sample maps exhibiting the matching circles to demonstrate that the whole computational chain (computation of the 
modes and implementation in a CMB code) is working. A complete and detailed study, along the same lines as the 
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FIG. 10: Contour plots of the function C{6i,92) for four realizations of the temperature anisotropies in a universe with a 
simply connected topology. When \9i — S2I is small, the correlation between circles exists. When \9i — ^2! is larger, the apparent 
correlation between circles comes from statistical fluctuations, which would be reduced for higher resolution maps. 



study done for the cubic torus in the previous section, will be presented in a follow-up article. 

As explained in Ref. because our universe is almost flat, observational methods such as the circles method will 
typically detect only a cyclic subgroup of the holonomy group, so the universe "looks like a lens space" no matter what 
its true topology is. It follows that lens spaces are particularly interesting to capture the observational properties of 
multi-connected spherical spaces. In particular, we showed l27l | that a cyclic factor Zp creates matching circles in the 
CMB only when — 1 > and that the second factor, if it exists, is in general undetectable. 

Let us emphasize some differences with the torus case. First, concerning the eigenmodes, let us take the example of 
a lens space L{p, 1) of order p. For p = 1 it reduces to and for p = 2 it reduces to projective space; more generally 
the index plays a role analogous to the size, L, of the torus in Euclidean space. The first non-zero eigenvalue is always 
J/ = 2 and has a multiplicity 3 for p > 2 and 9 otherwise. This constancy of the first eigenvalue contrasts sharply 
with the case of a cubic torus, for which the smallest eigenvalue scales as L~^. It can be understood by realizing that 
when p increases the space is becoming smaller only in one direction and remains large in perpendicular directions. 

The lens spaces L{p, 1) are globally homogeneous (like the torus) so that the coefficients do not depend on the 

observer's position (i.e., they are the same no matter where in the space you choose the bascpoint). Thus neither the 
correlation matrix nor the positions of the matching circles depend on the bascpoint. Unfortunately, this is not 

the case for a general lens space L(p, q). For a general lens space, the coefficients Cilm' correlation matrix C^^" 
and the positions of the matching circles all depend on the observer's position. For instance, the "canonical" choice 
of coordinates used in Sec. IIII B 21 for the toroidal coordinate system puts the preferred symmetry axes in the (xy) 
and {zw) directions (i.e., the axis are the intersection of S'^ with the (a;,y)-plane and the (z, w)-plane, respectively, in 
four-dimensional Euclidean space). From a cosmological point of view, this is a poor choice, because the observer's 
translated images are "atypically close". For example in i(12,5), which has cyclic factors Z3 and Z4, a generic 
observer will see three lines of four images each, but a nongeneric observer sitting on a symmetry axis will see a single 
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FIG. 11: Contour plots of the function C{6i,92) for four realizations of the Sachs- Wolfe part of the temperature anisotropics 
in a toroidal universe. No large correlations are found except when \9i — 62] is small and for the three pairs of matching circles. 



line of 12 images. 

For a globally homogeneous space L{p,l), the closest topological image is located at a distance x — 27r/p, so the 
topology is detectable just so p isn't too small. More precisely, the topology is potentially detectable if and only if 
T^/p < Xlss- For example, this implies that the topology is detectable for all p > 10 if fitot — 1 ~ 10~^. 

Circles match differently in a homogeneous lens space than in a tours. In a torus the circles match straight across 
because the holonomies are all pure translations. In a homogeneous lens space, by contrast, the holonomies are Clifford 
translations, and so the matching circles are still diametrically opposite but match with a twist that's a multiple of 
27r/p, because Clifford translations twist and translate of the same amount. 

Fig. 1151 shows a CMB map with resolution ^max — 30 for the lens space L(21, 1) considering the Sachs- Wolfe term 
only. A far more detailed discussion about CMB anisotropics in lens spaces will appear elsewhere 55] . 



VI. DISCUSSION AND CONCLUSIONS 



This articles describes the implementation of topology in CMB codes and gives explicitly the required tools to 
perform such an implementation in flat and spherical spaces. As emphasized in the Introduction, these two cases are 
observationally the most relevant for an almost flat universe. 

Examples of simulated maps were given in the two cases. Here we presented only low resolution maps due to 
the computational time limitation but higher resolution maps will be presented elsewhere. It was checked that the 
expected topological correlations (the matched circles) were present, confirming the quality of our simulations. 

Our method relies on the computation of the correlation matrix of the coefficients of the decomposition of the 
temperature fluctuation in spherical harmonics. This matrix encodes all the topological information. We emphasize 
that, due to the breakdown of global isotropy, this matrix is not purely diagonal. This also offers a working example 
to construct tests for the detection of deviation from global isotropy. 
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FIG. 12: Contour plots ol the function C{6\, 62) for four realizations of the Doppler part of the temperature anisotropics in a 
toroidal universe. Some correlation or anti correlation is found for two of the circles, but the signal is not very large compared 
to the statistical fluctuations. Note that the correlation would have been slightly larger if we would have simulated these maps 
from a higher resolution correlation matrix. 



We have illustrated the influence of different effects that will tend to blur these patterns and affect the perfect 
circle matching, namely the Doppler effect and the integrated Sachs- Wolfe effect. We also considered the effect of the 
thickness of the last scattering surface, but found it to be negligible on the scales considered here. A more detailed 
quantitative analysis of these effects on the detectability of the topological signal is left for future studies "sTl . 

A complete investigation of the detectability of the topology in coming CMB data requires the construction of 
reliable simulation tools. Besides the quantification of the amplitude of the effects cited above, one would also need 
to include all other observational effects such as instrumental noise, foreground contamination, etc. The present work 
paves the way to all these essential studies. 
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FIG. 13: Contour plots of the function C{6\,92) for four realizations of the integrated Sachs- Wolfe part of the temperature 
anisotropics in a toroidal universe. No significant correlation is seen on the matching circles, whereas large values are found for 
\C(9\,62)\ extending across broad regions as a consequence of the fact that the integrated Sachs- Wolfe effect appears on large 
angular scales. 



APPENDIX A: EIGENMODES OF CONSTANT CURVATURE THREE-DIMENSIONAL SPACES 



This appendix follows the work by Abbott and SchaefFer jsj and Harrison [Sg and borrows heavily from Appendix A 
of Ref. [l4f| . It summarizes, without proof, the explicit forms of the scalar harmonic functions solutions of the Hclmholtz 
equation |(2J). 

It is convenient to factor the eigenfunctions into radial and angular functions as 



y'fL{x,e,v)^Rfhx)Yr{o,v), 

with Yf^{9, (fi) being the spherical harmonics. The associated eigenvalues are K^ — k"^ — K, with 



With the normalization 



/ 



K <0 => ke[0,oo[ or ik e [0, V\K\], 

K = => ke [0,c»[, 

K>Q ^ k = {v+l)^/K , i^eN. 



where s/f (x) is defined in Eq. Q, the normalized radial functions take the form 

1/2 

'■kl 



ksKix) 



,-1/2-e 



-l/2+iu 



(cosh(/^x)) , 



(Al) 

(A2) 
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FIG. 14: Contour plots of the function C {61,92) for four realizations of the total temperature anisotropy in a toroidal universe 
(Sachs- Wolfe, Doppler, and ISW effects combined). The three pairs of matching circles at 61 = —62 — ±18°, ±39°, ±71° are 
still visible despite the presence of the Doppler and integrated Sachs- Wolfe contributions. 




with 

Lu^k/y/lKl, v = k/VK~l. (A9) 

In the case of spatially hype rbol ic spaces, this normalization is valid only for sub-curvature modes, and for the super- 
curvature modes {ik g [0, -y/jl?!]) the radial function is obtained by analytic continuation (see Ref. for details). 
The two numerical coefficients are given by 

e 

Mke = Y[ {{ly + If - n^) , = if £ > i^. (All) 

n=0 

For any function, we can perform the mode decomposition 

/(^) = E I ^5/2"^^- ^(^ + 1)2 } himyl^L ^ h^'n = j /(x):Vi?„/ ^^d^X, (A12) 
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FIG. 15: A realization of a CMB map in the case of the lens space L(18, 1). For purposes of illustration we chose a closed 
universe (fitot ~ 1-3) with a large last scattering surface (the radius Xlss is 0.88 times the curvature radius of the universal 
covering space S''), in order to produce higher resolution images with a large number of correlated circles. In each panel, the 
matching is obtained by performing a translation of nRc2n/18 and a twist of n27r/18 between the last scattering surface and 
its two copies. The six panels correspond to \n\ = 2, ... ,7, with each panel showing both the positive (n > 0) and negative 
(n < 0) translates. 



29 



choosing the sum or the integral according to whether the universal covering space is compact or not. The symbol 
stands for the square root of the determinant of the spatial metric. In the case of spatially hyperbolic spaces, 

the super-curvature modes add a term to this mode expansion, namely \k'^\d{ik) J2e m fkemykem^ 1^3 for 

details. 

In the spherical case, one can however find a solution of the Helmholtz equation l^l which does not involve Legendre 
functions. The radial part of the Helmholtz equation reduces, after setting y^f l^ — -Rfe^ (x)^/"(^7V')i to 



ix) dx 



d ^.2 d [S^ 

'■ke 



(fc2 -K)- 



i{e + i) 



ix) 



R'kV = 0- 



(A13) 



It is obviously much more convenient to work in a coordinate system where the curvature K reduces to 1. In terms 
of the dimensionless radial variable x defined in Eq. ^ , the Helmholtz equation then reduces to 



1 



sin X dx 



■ 2 - a n[S^ 



v(y + 2) - ^ „ _' 
sm X 



R 



kf. 



0. 



(A14) 



Note that this is a second order equation and that only one of the two independent solutions is well behaved at 
the origin, so the radial functions are completely determined once the normalization has been chosen. After setting 

-^i^^' ~ (sinx)^/i/^, it can be checked that it reduces to Eq. I|C11|) . the solution of which is simply given in terms 
of ultraspherical Gegenbauer polynomials as fy^ = Ai,iC^^^. The normalization condition (|A5|) implies, using the 
integral relation ljC12|) . that 



OTTITT)!' 



(A15) 



Expressing the spherical harmonics in terms of Gegenbauer polynomials by means of Eq. (|(]9(l , one ends up with an 
expression of the eigenmodes in terms of Gegenbauer polynomials only as 



sm' 



^l^?sin^X Cit\{cosx) 



COS! 



(A16) 



with 



Afm=Cm(2|m|-l)!!^ 



'(2^+1) {t~\m\)\ 



An {l+\m\)\ 

with C,ra given by Eq. IjGlOp and we used the notation (2|to| — 1)!! = (2|m| — l)(2|m| — 3) ... 1. 



(A17) 



APPENDIX B: CHANGE OF BASIS BETWEEN TOROIDAL AND SPHERICAL COORDINATES 

Sec, nil B 21 found the eigenmodes of lens and prism spaces in toroidal coordinates and converted them to spherical 

coordinates. In this appendix, we give the expression for the matrix gm^T mT necessary to perform the change of 
basis. 

The spherical coordinate system, as used in Eq. ((B)), is related to the embedding of the 3-sphere S'^ in four- 
dimensional Euclidean space by 

X — cosx, (Bl) 

y — sin X cos 0, (B2) 

z — sinxsin0cos</5, (B3) 

w = sin X sin sin (p, (B4) 

with 



< X < ^, (B5) 
< 6* < TT, (B6) 
< ^ < 27r. (B7) 
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The (complex) coefficients t ,„ mx characterizing the change of basis are defined by 

auimiTniT = +lf / Q,y£TmT(XT,6'T,</'T)3^i^^,'„*(X7^,'/')(sinx)^dx sin6'd6'd(p. (B8) 



In this expression the integer i ranges from to and m ranges from —I to while and range from —v to v. 
To compute this integral, one needs: 

1. to express the eigenmodes Q^i^. m-r a-s polynomials in x, y, z and w (see Ref. js^), and replace these rectangular 
coordinates by their expressions IBip . 

2. use the relations 

^^^|,,, {cos{\l^\e) ^ f 3?[(cosx + isinxcos0)l^-l] 1^>Q, 
\sin{\l.^\e) \ 9 [(cosx + ?sinxcos6')l^^l] < 0, 

sinl™-t XT [ = (sinxsinfl)'--! { '^o^d"^-)'^) ™- > (BIQ) 

I cos(|mT|iy9) I sin(|mT|(/3) < 0, 

3. develop the Jacobi polynomial appearing in Eq. (|57|l by using Eq. ICllll and with 

cos2xt - 1 = -2sin^Xsin^^. (Bll) 

This leads, after an easy integration on ip^ to the somewhat heavy expressions involving two sums [arising from the 
development of the Jacobi polynomials and the power in Eq. l|B9(l ] and an integral over 9 and x, 



and 



Here, the first and second line of the first brace are for > and < 0, respectively [see Eq. (|B9|l ]. the first and 
second line of the second brace are for ttit > and rriT < 0, respectively [see Eq. (jB10|l ]. and the numerical coefficient 
Cu t-T mT is given by 

_ r(|m^|+d+l) 

The function is explicitly given by 

.rdrn^l + 141 + 1) 



^0 V V 



(-1) 



11. 



r(|mT|+9 + l) 



X 



Y.^ ( ' j"' j ^ + I"^t|;^- |m|, |m| + 1 



2^ 

x/(|^t| -J,2g + .?+£+|mT| + (B14) 



Note that the index j is even when l^- > and odd when t^. < 0. This quantity involves only two sums, once the 
quantity /(p, q; n, a), defined by 

J(p, q- n, a) ^ x^{l - x^)'>/^C'::{x)Ax, (B15) 



is known. Using the expression (|C14p for the Gegenbauer polynomials in term of hypergeometric functions and the 
integral (|C16p . it can be shown that 



T( 9 ^ (-ir i3(g/2 + l,(p+l)/2) , ^ ^^^^ ^ ( ^ Hill ^\ 

I{p,q;2m,a) ^ —— — — — 1+ -If gi^a -m, m + a, — — ; -, ; 1 , B16 

2(a + m) B[a^7n + 1) \ 2 2 2 / 

7(p,,;2„, + l.o) . (_i).,. gW^+MP+-')/2) _ C_,„.„.^,^i.!L±j 3,i±£±i,i)(B17) 

ij(a, 771 + 1) \ 222/ 

where -B is the Euler Beta function. It follows directly from these expressions that j/(j/; ^, 777; 777^) = when 
V —\m\ + I^tI is odd. 
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APPENDIX C: SOME PROPERTIES OF SOME SPECIAL FUNCTIONS 

This appendix gathers some useful relations used in the article, to make the article more self-contained. 

The spherical harmonics IJ" are related to the associated Legendre polynomials P™ by (see Eq. (5.2.1) of Ref. 0|) 



vr^eM - y^l^ ^r(cos^) e--. (CI) 
They satisfy the conjugation relation (Eq. (5.4.1) of Ref. [58|) 

Yr{o, v) = (-l)"r,-™(^^, v) = yno, m 

the normalization (Eq. (5.6.1) of Ref. [s^) 

f sm0d6Yr*{0,'p)Y;r'(O,v)=Sei'5^„,,, (C3) 



/o Jo 

the closure relation (Eq. (5.2.2) of Ref. ^) 



e=o m=-e 

and the addition theorem (Eq. (5.17.2.9) of Ref. ^^I) 

^ Of -L- 1 

E Y,--{e,v)Yr{0',f') = ^^Pdcosa), (C5) 

m=-l 

where a is the angle between the two directions {9, ip) and {9' , ip'), and P£ is the Legendre polynomial. The Fourier 
transform of the spherical harmonics is given by (Eq. (5.9.2.6) of Ref. [SSf) 



r'dr d^ smed9j^-^Mk'r)Yri9,v)^\l^ r \, ^ ^(gk.^k), (C6) 
Jo Jo 



where ji is a spherical Bessel function, from which it follows that (Eq. (5.17.3.14) of Ref. 58]) 

oo 

g.k.r ^ ^(2^ + l)z^j,(fcr)P,(cos0k,r), (C7) 



t=0 

and (Eq. (5.17.4.18) of Ref. 



5^{v, - r,) . '.^^^1^ £ ^P,(eos^,,). (C8) 
The spherical harmonics can also be expressed in terms of Gegenbauer polynomials as (Eq. (5.2.6.39c) of Ref. js^ l 
Yr{9,^) = Cn.e-^^l^^l^^^i2\m\ - l)!!(sine)l™lcl':|:|/^(cos0) (C9) 

with C,„ defined by 

C„ = -J". (CIO) 

II m < 

The ultraspherical (or Gegenbauer) polynomials C" are solutions of the differential equation (Eq. (22.6.5) of Ref. 59]) 

(1 - x^)y" - {2a + l)y' + n{n + 2a)y = 0, (Cll) 
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and they satisfy the normalization condition (Eq. (7.313) of Ref. |60l |) 

L n\(n + a)[r(a)\ 

if 3?(a) > -1/2. 

The Jacobi polynomials are given by (Eq. (22.3.2) of Ref. [s^) 

^" n!r(a + ;9 + n + l)^^l^jJ 2^ T(a + j + 1) ' ^^^^^ 

under the conditions a > — 1 and /3 > — 1. Interestingly, the Gegenbauer polynomials can be expressed in terms of 
hypergeomctric fmictions as (Eqs. (8.932.2,8.932.3) of Ref. 0|) 

(—1)" 

CLi^) = (, + 4^(1,, + !) F i-n, A + n; 1/2; , (C14) 
^2"+!^"^^ = -g^^=^^(2a;)F(-n,A + n+l;3/2;x2), (C15) 



which satisfies the integral property (Eqs. (7.513) of Ref. |60j ) 

x'*-i(l - x'^yP {-n, a- b; x^) dx = ^B{v + 1, s/2) 3F2 (-n, a, s/2; 6, + 1 + s/2; 1) (C16) 
if 3fi(s) > and ^{v) > -1. 



MAP homepage: 'http : //map . gsf c . nasa . gov/ 1 . 



Planck homepage: http : / /astro . estec . es a.nl/Plaack/] . 
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